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I. INTRODUCTION 



It has been known for many years [jjj that next-to-leading order (NLO) calculations in 
perturbative quantum chromodynamics (pQCD) for cross sections of various hard scattering 
processes, lead to functions singular at the edge of the relevant phase space. For example, for 
Drell-Yan (DY) process we define z = Q 2 / s where Q 2 is the invariant mass of the produced 
lepton pair squared and s is the total partonic invariant mass squared. At nth order in 
pQCD we encounter terms like — z)^ 1 ln 2n_1 (l — z)\. In the limit z — > 1 these singular 
functions — when Mellin transformed — give rise to large double logarithms in moment space 
(a s ln 2 iV) n . When N is large, the product of these large logarithms with a small coupling 
constant a s is not necessarily small and an all-order resummation in perturbation theory 
is needed to make predictions reliable. The purpose of a resummation is to control these 
large logarithms and to obtain a closed expression that sums part or all of these logarithms. 
For DY and deep-inelastic scattering (DIS) processes, an extensive theoretical work was 
performed in Refs.0, HI, according to which, one has to establish first a factorized form of 
the cross section into a hard part ("Wilson coefficient"), jet part and soft part, and then 
make a clever use of an evolution in parton rapidity. We will have more to say about this 
in section IV. 

Recently an effective field theory approach was implemented [H to re-sum leading logs 
(LL) and next-to-leading logs (NLL) for DIS at large momentum transfer and near the 
threshold limit x — > 1 (where x is the longitudinal hadron momentum fraction carried by 
the quark). That analysis was performed within the soft-collinear effective theory (SCET) 
framework 0,0,0]. By utilizing this effective theory for the inclusive DIS (for large momen- 
tum transfer and near threshold) cross sections, one can re-derives the full QCD results (as 
a consistency check of the effective theory!) More importantly, the exponentiation of the LL 
and NLL in moment space appears quite naturally in this context following the solution of 
a simple RG equation 0]. 

In this work, we extend the application of SCET to the DY cross section in the threshold 
limit z — > 1 and show how the resummation is obtained. We work in the center-of-mass 
(c.o.m.) frame of the incoming quark and anti-quark, and identify three relevant scales: Q 2 , 
Q 2 (l — z) 2 and Aq CD . In the threshold limit we have 1 - z « 1. Assuming Aq CD /Q 2 <C 
(1 — z) 2 , we have a clear separation of scales. Given these three scales one has to perform 
a two stage matching and running: First, at Q 2 we match the full QCD (annihilation) 
current onto the SCET current, giving rise the so-called SCETj in which virtuality of order 
Q 2 is integrated out. Then we match, at the intermediate scale Q 2 (l — z) 2 , the SCET DY 
cross section onto a product of SCET parton distribution functions (PDF) 8], giving rise 
to SCETn in which virtuality or order Q 2 (l — z) 2 is integrated out. Here we assume that 
the usoft interactions cancel as in the standard Drell-Yan factorization which is valid in 
the kinematic region of our interest 0. The second matching yields the infrared-safe hard 
scattering coefficient function. Then, the PDF will run down in scale and be matched onto 
local twist-two operators. 

This paper is organized as follows. In section II we give the basic notation, definitions and 
power counting of SCET. In section III we perform the two-stage matching and write down 
the anomalous dimensions. We compare our results for DY (as well as those of DIS) with 
the full QCD calculations. In section IV we derive the resummation (and exponentiation) 
of the large logs and compare our results with those given in 0, HI . We conclude the paper 
in Section V. 
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II. SCET PRELIMINARIES 



Soft-collinear effective field theory was invented to handle processes with multiple mo- 
mentum scales 0, 0, Q • This is the case in DY and DIS when final-state hadrons have an 
invariant mass much greater than Aqcd but much smaller than the highest momentum scale 
of the problem. In this section, we briefly describe some of the concepts and notations in 
SCET. We choose 

n" = -L(l, 0,0,-1), n" = -±=(1,0,0,1), n 2 = n 2 = . (1) 

For a generic four-vector I let = ^(^° ± I 3 ) and write / = (Z + , l~ , l±). Then 

1+ = n -l, r =n-l, l 2 = 2l + T + l\ = 2n ■ In ■ I + l\ . (2) 

With this notation one has to trivially modify the SCET Feynman rules given in 0|. 

SCET describes the interaction of "collinear" quarks with "collinear" and/or "soft" glu- 
ons. If we have a fast-moving quark or gluon in the n-direction then we assign to its 
momentum components the following scaling: p = (p + ,p~ ,p±) ~ <3(A 2 , 1,A) where Q is 
some relevant hard scale and A is a small parameter identified by considering the different 
momentum scales available. A quark with such an assignment is called n-collinear and is 
described by £ n field and n-collinear gluon has field. Similarly, for n-collinear quark field 
we assign p = {p + ,p~ ,p±) ~ Q(l, A 2 , A) and is described by ^ field. For the DY case, and 
in the c.o.m. frame, we let the incoming quark be moving in the +z direction (i.e Pi is 
"large") and the incoming anti-quark moving in the —z (i.e. P2 is "large"); then these fields 
are described by £ n and £ n respectively. At the scale Q, after integrating out virtuality of 
order Q 2 , the electromagnetic current in SCETi is given by 

f = C{Q) lnW n ^Wk n . (3) 



The coefficient C (Q) has been calculated in and will be presented in the next section. The 
W n and Wl are the familiar path-ordered Wilson lines and are required to insure collinear 
gauge invariance of the current operator This Wilson line is identical to the "eikonal" 

line introduced in 



W n (x) = P exp 



X 



ig ds n ■ A n (sn) 



(4) 



where the covariant derivative is: = — igT a A a ^. A "soft" gluon has the momentum 
scale: k = Q(X 2 , A 2 , A 2 ). A SCET PDF for a n-collinear field is defined as the matrix element 
of the operator 

1 f°° r i 

O q (r + ) = — J dze~ lzr+ [lnW n \nz) ft |w&nj(0) , (5) 

between n-collinear quark fields. It is clear that for DY process one needs to consider the 
n-collinear version of Eq. (5), i.e., interchanging n n and taking the matrix element 
between n-collinear fields. 

In this work we use the Feynman gauge and the MS scheme in d = 4 — 2e. We will also 
use on-shell dimensional regularization (DR) to regulate both the ultraviolet (UV) and the 
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infrared (IR) divergences. The choice of IR regulator in SCET is a delicate matter and 
one has to make sure that the IR regulator accurately generates the IR behavior of the full 
theory. In Ref. quark off-shellness was used as IR regulator, however all the results can 
also be obtained by using pure on-shell DR. This was checked explicitly. We will also follow 
the normalization conventions for the DY and DIS cross sections given in []] . Thus the Born 
cross sections are 6(1 — z) and 6(1 — x) respectively. 



III. MATCHINGS AT Q 2 AND Q 2 (l - z) 2 FOR DY: pp -» l+l~X 

In this section, we derive the matching conditions for DY at Q 2 and Q 2 (l — z) 2 at one- loop 
order. From these, we re-derive the known one-loop coefficient function in SCET. To make 
the discussion self-contained, we first review the result for DIS process 0. 



A. Review of DIS in SCET 



Since the DY and DIS processes are related to each other, let us review, briefly, how the 
analysis for DIS in SCET was performed. Choosing the Breit frame, let the incoming quark 
be n-collinear, (i.e. moving in the +z with large pf and small p\) and the outgoing quark 
be n-collinear, and let x = Q 2 /2pi ■ q. The invariant mass of the final hadronic state for DIS 
in the limit x — > 1 is: p 2 x = Q 2 (l — x). Simple kinematic considerations show that for either 
collinear or soft gluon radiated into the final state: p\ = Q 2 \ 2 + 0(X 4 ) thus we identify: 
A 2 ~ 1 



x. 



Io dx 



x 



N 



l f(x). M N 
oo, thus one 



The Mellin transform of a function / is defined as: M^[f(x)\ /(| 
is also known as "N th -moment" of /. The limit x — > 1 corresponds to N 
naturally, identifies, 1 — x ~ 1/N. 

At scale Q 2 , the final hadronic states may be effectively considered as a massless jet 
moving in the n direction. One matches the full QCD current (i.e. only incoming and 
outgoing quarks) onto SCETi current given in Eq.(3). The matrix element of the full QCD 
current in on-shell DR is given by 



(P2\f\pi) = *f\l + -^C F 



O's 



'^-^l 21n ^ + 3 J- ln ^- 31n ^- 8 + T. 



6) 



The 0(a s ) virtual corrections to SCET current are scaleless and vanish in pure DR, thus 
the matching condition at Q 2 equals the finite part of Eq. (6) 



Cdis(aO 



1 + > 



, 9 /i 2 , U 2 7T 2 

- ln ^- 31n ^- 8 + T 



(7) 



The anomalous dimension of SCETi current satisfies 

dC D i S (u) 



du 



7i(/i) Cdis(«) 



which gives 



(8) 



7i(M) = ~^C F 



41n^ + 6 



(9) 
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As the intermediate scale Q 2 (l — x), the final-state invariant mass p\ must be considered 
as "large", and the jet-like final-state n-collinear hadronic modes can be integrated out. 
This is done by matching a product of SCETi currents onto an effective theory (SCETn) 
where the n-collinear quark is "replaced" by a soft Wilson line, thus one is led to the matrix 
element of SCETn quark operator given in Eq. (5) between quark states. This is very much 
like the standard operator product expansion. If one computes the matching condition at 
Q 2 (l — x) by using pure on-shell DR, then one has to consider only the contribution of 
Feynman diagrams where a real gluon is emitted into the final state, because the one-loop 
virtual diagrams are scaleless and thus vanish. 

If one computes the real-gluon emission diagrams in SCETi, one finds the 0(a s ) matching 
condition (coefficient function) in 

M ms (x) = ^C F d(0<x<l)\2 

^f-^f + i-yK-*)}- < 10 > 

The above result equals exactly to the finite part of the full QCD calculation of the real 
gluon emission diagrams 0, up to terms that vanish in the large moment limit. Taking the 
product, Cj) IS (Q)[8(l — x) + M.ms( x )] yields the full QCD coefficient function. Therefore, 
the usual one-step matching in DIS is divided into two steps, virtual and real corrections, 
in SCET. This is also true in Drell-Yan process to be discussed below. 



ln(l — x) 


+ 




(1-z) . 


4- 


L fi 2 21 



B. Drell-Yan Process 

In the Drell-Yan process, the invariant mass of the final hadronic state is 

p\ = Q 2 (l + i - 1 - 1) s Q 2 [(1 - Xl )(l - x 2 )} - Q 2 A 4 , (11) 

\ Z X\ X2 J 

where 

Q 2 Q 2 Q 2 , 19 , 

Z = , X\ = , X2 = • (12) 

s 2p x - q 2p 2 ■ q 

Since 1 — z ~ A 2 , p\ ~ Q 2 (l — z) 2 . To integrate out the final hadronic states, we match at 
the intermediate scale Q 2 (l — z) 2 . 

At the scale Q 2 , the final-hadron states effectively have p 2 x = 0. We match the full QCD 
annihilation current: qq — > 7* onto the SCETi current. The matching condition can be 
obtained simply by analytically continuing the full-QCD matrix element in Eq. (6) from 
space-like q 2 = —Q 2 to time-like q 2 = Q 2 , 



a 



•4— L( 2 ]n£ + 3 
em em V Q 2 



(13) 
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The matching condition at Q 2 is thus 

a 



1 + 



-In 



Q 2 



, /i 2 7tt 2 



(14) 



The anomalous dimension is the same as Eq. (9). It should be noted that eiR-dependent 
terms in Eq. (13) equal the negative of the euv terms in the effective theory. This is so 
since, as mentioned earlier, the SCET virtual contributions vanish in pure on-shell DR due 

to cancellation of the form — — , and because the IR behavior of the effective theory 

must match that of the full QC$D. Tnus the UV counter-term for the current in the effective 
theory equals 



c.t. 



a, 
27r' 



C f 



_ 1 _ 3 1 i n V* 



2e 



Q 2 



(15) 



Below Q 2 , the DY process is described by SCETi. Because the virtual diagrams vanish 
in on-shell DR, we need only consider real emission diagrams given in Figs. 1 and 2, where 
the cut lines go through the real gluon in the final state. Let us take, in the c.o.m. frame, 
the momentum fc M of the emitted gluon as: 

(16) 



fc" = (|fc|,...|fc|cos0) 
and introduce the Mandelstam variables (besides s) 



t 



- — (1-*)(1 

z 



y), 



u 



-£<1 

z 



z)y 



(17) 



where: y = |(1 + cos 9). The calculation of the contributions from the above diagrams 
is performed by taking the square of the amplitude for emitting a real gluon and then 
integrating over the phase space for the production of a massive photon p| 



/4tt 
8tt VQ 2 



PS = — — 



vl-2e 



r(i 



dy(y(l -y))' 



Diagrams (c) in Fig.l and (a)-(b) in Fig. 2 are zero due to n 2 



n 



(18) 

0. Diagrams (a) 

and (b) have the same contribution as their full QCD counterparts. The sum of the three 
remaining contributions completes the total result to that of full QCD which can be found 
in, e.g, jl(. We perform the phase space integration for the sake of completeness, keeping 
only singular terms in the limit z — > 1, 



1by(z) 



^M(0<-€ 1V/(1 " Z 



r(i 



ir/ w -»)rx 1 1 



a 



—C F 9(0 < z < 1) x 2 
2tt 



-7?*(1-*) 



1 - Z 



+ 4 



ln(l - z) 
1 - z 



. U 2 n / 1 , o U 2 7T" 

1 + £ l n ^_ + £ ^_ b ^-_ 



a 



—C F 9(0<z< 1) 
2n 



2 2 u 2 

- 5{l-z) + -\n^-5{l-z) 



Q 2 



e 1-2 



1 



2 ^ 



7T" 



(19) 
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in % 



in 



(g\/\/V 



'\ZV\igl 



(a) 



in 



(b) 



(d) 



(e) 



FIG. 1: Collinear real-gluon contributions to DY crosssection in SCET. Diagram (c) represents 
the two possibilities of n- collinear and fi-collinear gluons. Diagrams (d) and (e) have mirror image 
graphs. 



in n 



,®\/\/V 



in '' 



in N 



in '* 



(b) 




FIG. 2: i?eaZ soft-gluon contributions to DY cross section in SCET. Diagram (c) has mirror image 
graph. 



In the above calculation we ignored the contribution from lnz/(l — z) since this function 
is regular in the limit z — > 1 and the corresponding Mellin moments vanish at large N. As 
expected, the result is identical to the full QCD calculation of the real gluon contributions 
to the DY differential cross section in the limit z -> 1 0. The e-terms in Eq. (19) are a 
combination of UV and IR contributions. 

To see how to interpret these terms, we need to consider the PDF given in Eq. (5). In 
pure on-shell DR, the 0(a s ) corrections to the unrenormalized parton distribution vanish 
due to scaleless integrals. The UV-behavior of the SCET PDF is given by 0]: 2 " e s P q ^- q (z) 
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where in the limit z 



P, 



q^q 



'1-z) 



0(0 <z < 1) , 



(20) 



which is the well-known Altarelli-Parisi splitting kernel. Thus the renormalized parton 
distribution is — -^^Pq^q( z ) ■ Writing the e-terms in Eq. (19) as 



a 



-^C F 9(0<z< 1) 

Z7T 



2 3 2 u 3 , . , 



1 - z 



-6(1 - z) 



Taking into account the counter-term for the SCET current, Eq. (15), and matching the DY 
cross section onto a product of two PDFs, all the e-dependent terms cancel. 

The finite part of Eq. (fT9")) yields 0(a s ) matching condition (coefficient function) 



ft.. 



M DY (z) = ^C F 6(0<z< 1) 



7T 

2 ~ ~2 



/ ln(l - z) 

\T^z~ 

2* 



1 - Z 



(21) 



The moments of the matching condition at Q 2 (l — z) 2 , Eq. (21), are 



M D y(N) 



a i 
2n' 



C f 



In 2 



~Q 2 ~ + ~6 



(22) 



where N = N exp(7^) and 7^ is the Euler constant. In order to minimize the large logs we 
choose fi = Q/N. This is different from the DIS case where at the intermediate scale we set 

H = Q/vN. Thus for DY we identify 1 — z ~ 1/N. At this stage it is appropriate to clarify 
the following point: the above calculation demonstrates the all order factorization theorem 
°f ^e DY process to first order in a s , i.e, all the IR divergences are contained into a 
product of two PDFs defined in Eq.(5). 



Below the scale Q 2 (l 



we have SCETn in which the all hard scales have been 



integrated out. The running of the moments of the PDF with /i is governed by the well- 
known anomalous dimension —72,7V [3 03 where 



72,iV 



2tt 



N 1 

4 Et 



j=2 



3 N(N + 1) 



(23) 



In 0] it was shown, to first order in a s , that the running of the PDF comes from Fenyman 
diagrams where only collinear quarks and gluons interact and the soft contribution (or usoft- 
in the terminology of Ref.0) vanishes. This is of-course consistent with the factorization 
theorem of DY which we assume that it holds to all orders in perturbation theory. 
In the large N limit we have 



a 



72.JV 



2vr 1 J 



(24) 



At the lowest scale un, one identifies the PDF moments as the matrix elements of local twist 
operators A N (fj, ) |l4| . 
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IV. RESUMMATION AND EXPONENTIATION 



By now, we have obtained all the ingredients we need to write down an expression for 
the moments of the DY differential cross section (^jr)- By taking the square of the current 
matching coefficient at fi 2 = Q 2 , running down the scale with anomalous dimension 271 to 
/i 2 = Q 2 /N , we then multiply with the moments of matching condition into the product of 
the matrix element of PDF. Then we run down the scale to n ~ Aqcd with 272 and match 
onto the moments of local twist-two operators. In summary, one gets 



da 



DY 



dQ 



) = C 2 Y (QK /DY1 (1 + M DY (N, l i = Q/N)) e-^A 2 M , (25) 

/ N 



where 



( DY1 



f Q dfi , s [ Q / N dn , N 

/ __^2 7 i(/i); / DY2 = / — 272M, (26) 



-Q/N 
' Mo 



where 71 (« s , Q 2 /fi 2 ) is the anomalous dimension of the SCET current and 72(0^, N) is the 
anomalous dimension of the twist-two operator. 
The above can be compared with the DIS result, 

F 2N (Q 2 ) = C 2 IS (Q)e- /DIS1 (l + M ms (N, fi = Q/t/N)) e'^A^o) (27) 

where 

271 (/i); J D is2 = / 
Iq/Vn V J no 



'DISl 



/ _-^2 7l ( A i); /dis2=/ — 7 2 (/i). (28) 

J Q/v N V Jno V 1 



Let us define a physical observable which is a ratio of the moments between DY and DIS 
squared, 



(da DY /dQ 2 ) 



A N = v — 7" > N . (29) 

-To 



2N 



At one-loop order, one has 

A N = l + —C F (2\n 2 N-3\nN+l + -n 2 ), (30) 
2n \ 3 / 

which is consistent with the known result. When summing over higher-order corrections of 
form a™ ln m N with m < 2n, Ajv can be expressed in terms of an exponential form, exp(/), 
with / = ln]V/_i(a s ln]V) + f (a s \nN) + a a /i(a a \nN) + .... The first term yields the 
leading double-logarithm resummation, and the second gives the next-to-leading logarithm 
resummation (NLL), etc. In the present calculation, we can have an accuracy up to NLL. 
To NLL resummation, A^v can be expressed as 

r Q da r Ql ^ da 

\nA N = 2ll (a s (ri)-+l 2[ 7l (a s (/i))-72(« s (A0)]-, (31) 

Jq/vn H Jq/n ^ 



9 



where the anomalous dimensions have the form 



7 2 = 2r 1 lnA7 + r 3 , 



(32) 



to all orders in perturbation theory. To NLO, one has T 3 = 2r 2 , and the above expression 
can be expressed as 



In A 



N 



2 1 1 it 

In- 1 y 



^ Q da 
UyQ V 



(33) 



This result does not have an infrared renormalon or Landau singularity because the expres- 
sion is manifestly valid only when the scale Q/N is much larger than Aqcd- If one takes 
the limit N —>■ oo with a fixed Q, the above result can also be expressed as 



An = exp 




dk 



(l-z) 



T 1 (a s (k 2 ))— + T 2 (a s ((l - z)Q 2 )) 



(34) 



which is the same as what has been obtained by Catani and Trentadue |3j , and Sterman 
This expression has the infrared renormalon problem because the integral covers the soft 
momentum region where the strong coupling constant is ill-defined. 

It is instructive to recall how the resummation in the large- x limit is achieved in Sterman's 
approach 0. The DY cross section and DIS structure function can be factorized in x — > 1 
limit into a product of hard parts, a soft part, and collinear parts. Different parts may be 
gauge-dependent, but the product is not. The renormalization group equations for individual 
parts can be derived straightforwardly. 

The hard part contains lines with momentum of order Q and is included in the first stage 
matching in SCET approach. The soft part is defined as the matrix element of Wilson lines 
and is infared finite (for a definition in SCET, see 0]). This soft contribution comes in fact 
from gluon lines of virtuality of order Q 2 (l — x), which can be calculated in perturbation the- 
ory when Q 2 {1 — x) ^> Aq CD 0]. In SCET approach, this contribution is taken into account 
by the second stage matching. The collinear parts are either special parton distributions or 
jet functions or both. They depend on a rapidity cut-off and obey an evolution equation in 
rapidity. This equation is quite general, and when solved, contains the resummed double 
logarithms. When Q 2 (l — x) ^> Aq CD , the special parton distriubtions can be factorized in 
terms of ordinary parton distributions, whereas the jet functions can be calculated entirely 
inpQCD. 



V. CONCLUSION 

Following a previous work on summing large logarithms as x — > 1 in DIS, we use the 
same SCET formalism to study the resummation in DY. The steps are quite similar and the 
result is shown in Eq. (|25|). The result is consistent with the known result, but is free of the 
Landau singularity. 
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